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Abstract 

How to effectively solve the eigen solutions of the nonlinear spinor field equa- 
tion coupling with some other interaction fields is important to understand the 
behavior of the elementary particles. In this paper, we derive a simplified form 
of the eigen equation of the nonlinear spinor, and then propose a scheme to 
solve their numerical solutions. This simplified equation has elegant and neat 
structure, which is more convenient for both theoretical analysis and numerical 
computation. 
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1 Introduction 

Almost all the elementary fermions have spin-|, which can be naturally described by 
spinors, so today spinors and spinor representations play a more and more important 
role in mathematical and theoretical physics. Noticing the limitations of the linear field 
equation, many physicists such as H. Weyl, W. Heisenberg, once proposed the nonlinear 
spinor equations[l, 2, 3, 4, 5, 6] to construct a unified field theory for elementary 
particles. However they have not gotten many definite results due to the mathematical 
difficulties. The rigorous solutions for some simple dark nonlinear spinor models were 
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obtained in [7, 8, 9, 10], and we found it can provide negative pressure to guarantee a 
singularity-free and accelerating expanding universe [11, 12]. 

The theoretical proof about the existence of solitons was investigated in [13, 14, 
15, 16, 17, 18]. The symmetries and many beautiful conditional exact solutions of the 
nonlinear spinor, vector and scalar differential equations are collected in [19] . However 
lots of these exact solutions seem to be non-physical. 

The spinor with its own electromagnetic potential was studied in [20, 21, 22, 23, 24], 
it was disclosed that the nonlinear spinor equations have particle like solution with 
anomalous magneton, and imply the exact classical mechanics and quantum mechanics 
for many-body[25, 26, 27]. 

In this paper we derive a simplified form of the eigen equation with general meaning 
for nonlinear Dirac equation, and then give a scheme to solve the solution. Denote the 
Minkowski metric by rj^ = diag[l, —1,-1,-1], Pauli matrices by 

: :)} 

In this paper, we adopt the Hermitian matrices (1.2) instead of Dirac matrices 7 M , 
because this form is more convenient for calculation. 

For the system of a nonlinear spinor field <f> in the 4 — d potential A 11 , the Lagrangian 
describing the motion is generally given by 

C = + K(^<9 M - eAp) - /ic 7 ]0 + F(7, /3), (1.3) 

where /i > is a constant mass, F are the nonlinear coupling potential, which is usually 
the even polynomial of 7, $, and 7, $ are the quadratic scalars of defined by 

7 = + 70, $ = 4> + P<t>- (1-4) 

We can check 7 is a true-scalar, but f3 a pseudo-scalar. 

The variation of (1.3) with respect to gives the dynamic equation 

a^Hid^ - eA^ = (^c 7 - F l7 - F^)0, (1.5) 

where F 7 = Fp — |£. In the Hamiltonian form we have 

hid t (f) = H(f), H = a- {-hiV - el) + eA + (/ic - F 7 ) 7 - F p p. (1.6) 



[a 3 : 





Define 4x4 Hermitian matrices as follows 
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In this paper we denote A = (A 1 , A 2 , A 3 ) to be the spatial part of a contravariant 
vector N 1 . 

It is easy to check that the current conservation law holds = 0, so we can take 
the normalizing condition as follows 



\ 2 d 3 x=l. (1.7) 

R 3 

The 4 — d potential produced by spinor itself takes the following form 

d a d a A» = ea» = e<$> + ot<S>. (1.8) 

2 Simplification of the Equation 

Consider a spinor keeping motionless in an external magnetic field B = (0, 0, B). since 
the scale of the elementary particle is very small, we take external field B as a constant. 
By V x A ext = B, we have the general solution for external vector potential 

A ext = ^B(-y,x,0) + V$, (2.1) 

where <&(x) is any given smooth function. 

In the spherical coordinate system (r, 9, (p), we have 

a ■ V = a r d r + (aed e + a^d^), (2.2) 

where (oy, ag, a^) is given by 

cos6> sin0e-^\ if-smO cosfle - ^ 1 / -ie~^\ 1 
v sin Be* - cos 9 J'r \cos 9e^ sin 9 J'r sin 9 \i e ^ / J ' 

Let J be the angular momentum operator for the spinor field 

j = f x (-MV) + -hS, S = diag(<r,a), (2.4) 

then any eigenfunction of J 3 = —hid^ + ^hS 3 takes the following form 

K^i — til (2.5) 

with (« = 0, ±1, ±2, ■ • • ), where Uk, Vkik = 1,2) are functions of r, 9 but independent 
on (p and t. In this paper the index T stands for transposed matrix. 

For any spin-| particle, it has a pole axis. If we set the pole axis as coordinate 
x 3 = z, then J 3 is commutative with the nonlinear Hamilton operator (1.6) by a £7(1) 
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gauge transformation for spinor as e* l (/>, which removes the uncertain function from 
external vector potential (2.1), thereby we have 

A ext = ^B(-y, x, 0) = i.Br sin 0(- sin <p, cos ip, 0). (2.6) 

For the above symmetric form of vector potential (2.6), substituting (2.5) into (1.6) we 
can check that all functions Uk,Vk can take real number. For the vector potential pro- 
duced by (ft itself, we will find below it also takes the form of (2.6). This simplification 
of <ft may be the essence of the gauge symmetry. 

In what follows, we set h = c = 1 as units for convenience. Making variable 
transformation 

u = ui(r,6) + u 2 {r,6)i, v — vi(r, 9) + v 2 (r, Q)i, (2.7) 

then we have 

d = \u\ 2 + Id 2 d = (uv — uv)i, 

11 _ _ _ ' ' (2.8) 
7 = |w| 2 — \v | 2 , $ = uv + uv, 

with a = d(— sin <p, cos tp, 0). Substituting (2.7) and (2.8) into (1.3) we get the La- 
grangian of the eigen states as follows 

C = Re[e di (u(d r + i r d e )v-v(d r + i;de)u)] - _^(« + !)(««- ««) 
+ {m — eA )(| , u| 2 + \v\ 2 ) — ie(uv — uv)A — fi{\u\ 2 — \v\ 2 ) + F, 

where 

A = -(fx A) ■ e z = (cos ip A y - siny?^) = A(r, sin 6) (2.10) 

including both external and inner vector potential. By variation with respect to u, v, 
we get an elegant equation with double-helix structure 

e ei {d r + l -d e )u = ^[(k + \)u -\u} + {u + m- eA - F^v + F?u + teAu, 
e ei (d r + id e )v = -A- 9 [( K + \) v -\v\ + (iJL-m + eA - F 7 )tx - Fpv + ieAv. 

By (2.11) we find that, k — corresponds to spin | and k = — 1 corresponds to spin 
— |. Generally the coordinates r and can not be separable for nonlinear equation. 
The energy functional for (2.11) is given by 

E = 2n J r 2 sin 9drd8 {—Re [e 8i (u(d r + ^d e )v — v(d r + ^d e )u)~\ — F 

+ bb( K + I) + eA ] Kuv - uv) + eA {\u\ 2 + \v\ 2 ) + fi{\u\ 2 - \v\ 2 )} . 

The eigen solution of (2.11) is just the extreme point of E under the constraint of 
normalizing condition 

2tt J{\u\ 2 + \v\ 2 )r 2 sin9 drdd = 1. (2.13) 
(2.13) is also the quantizing condition of the energy spectrum[9, 10]. 
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3 A Scheme for Solving Solution 



For general potential and F, the analytic solution of (2.11) u and v can not be 
solved. However they can be conveniently expressed by Fourier series of 6, and the 
equations of the radial functions can be derived by variation principle. For any given 
integer N > 0, define 2N + 1 vectors 

r(0) = (e-^e- 2 ^- 1 )*... ,e 2 ^ ei ,e 2mi ), (3.1) 
U(r) = (?7_ JV (r),i7_( JV _i)(r),...,^ (JV _i)(r),i7 JV (r)) T , (3.2) 
V(r) = (y_ JV (r),y_ (JV _ 1) (r),...,V( JV _ 1) (r),Vi v (r)) r . (3.3) 

The eigen solution of (2.11) with even parity must take the form 

N N 

u = T-U= U n e 2ne \ v = T- Ve ei = ^ V n e^ 2n+ ^ 6i , (3.4) 

71= -JV 71= -JV 

and the eigen solution with odd parity takes 

N N 

U = T- Ue ei = U n e {2n+1)ei , v = T-V = J2 Ke _2nW . (3.5) 

n=-N -N 

In what follows we only consider (3.4). For (3.5) we have similar results. 

For the cases k^O and k ^ — 1, i.e. for the cases with nonzero magnetic quantum 
number, the solution must have consistent conditions at 6 = 0, n as 

JV JV 

u(r, 0) = u(r, n)= U n( r ) = °> <r, °) = n ) = K ^ = °" ( 3 - 6 ) 

n=-N n=—N 

For this case, (3.4) minus (3.6) we get 

JV TV 

u = ^ [/ n (e 2 ^-l), v = V n {e~ 2n6i -l)e ei . (3.7) 

n=-N n= -JV 

By the form (3.7) we have 

= (i + e 2« + e m + . . . + e 2(»-i)«) e « (3.8) 

which removes the singularity of (2.11) at = 0, 7r. 

For the spin | state, i.e. for the case k = 0, we can check from (2.11) that the 
solution U, V are all real functions. But for the spin — \ state, i.e. for the case k — — 1, 
the solution U, V are all pure imaginary functions. In what follows we only consider 
the real case. For the covariant quadratic forms (2.8), by (3.3) we have 

a = U T PU + V T PV, a = U T (Q + -Q)Vi, 
7 = U T PU - V T PV, (3 = U T (Q + + Q)V, 



7 = 



a = 



where P = T + T and Q = r T re e * are (2N + 1) x (2iV+ 1) matrices with components as 

P m , n = exp[2(ra-m)0i], (-N < n,m < N), 
Q m;n = exp[(2(n + m) + l)6i\. 

By (3.9) and (3.10) we have 

N 1 

*° = E 2^^ m + Kym ^ e " 2(n " m)9i + e2(ri " m)ei )' ( 3 - n ) 

n,m= — N 

^ 1 

x; 2 (Kt ^" Kym)(e " 2(n " m)w+e2(n " m) ^ ) ' (3 - i2) 

n,m=—N 
N 

U nVm(e~ 2(n ~ m)ei ~ e 2{ - n - m -^ ei )e e \ (3.13) 

n,m=—N 
N 

= U n V m {e- 2(n - m)ei + e 2(»-™-i)W) e <M. (3. 14 ) 

n,m= — N 

The dynamic equation of the potential A is given by 

— AA = ea = e[a (r) sin 9 + a\{r) sin 30 H ], 

= e[a (r)(e- m -l) + a 1 (r)(e- 4ei -e 2f?i ) + ---]e^. (3.15) 

Substituting (3.4), (3.7), (3.8) and (3.11)~(3.15) into (2.11) and directly comparing 
the coefficients of all e 2ndl , we can easily get a truncated ordinary differential equation 
of U(r),V(r). However the convergence this cut-off equation to the original solution 
needs proof. A more credible method to get the efficient equation of U(r),V(r) is 
via variational principle. The variational equation can be obtained by the following 
procedure. Define operators 

f u = \ f T T (6)e- dl sin 9d9, T v = \f T + (9)sm9d9, (3.16) 

then T u left multiplying the first equation of (2.11) and T v left multiplying the second 
give the variational equation of U(r), V(r). The coefficient matrix of U'(r) and V'(r) 
is the same positive definite symmetric matrix with components 

i r 1 

M nm = - / P nm sm9d9 = ; -. (3.17) 

n,m 2 ^ n,m X _ 4 ( n _ m )2 V ) 

The other coefficient matrices can also be similarly obtained. 

The convergence of expansion (3.4) and the consistent condition (3.6) seem to 
have closely relation with the structure of nonlinear potential P(7, $), which reflects 
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the important properties of the elementary particles such as the exclusion princi- 
ple. The above procedure is valid for extensive models. The neat and elegant re- 
sults are profoundly rooted in the quaternionic structure of the physical variables and 
spacetime[28, 29, 30], so the 3 + 1 dimensional Universe is a miraculous masterpiece 
with unique feature. 
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